Analysis of Hadronic Properties at SPS energies from a Statistical Model 
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Using a statistical model of multiparticle production by Chou, Yang and Yen, which looks similar 
to the thermal model and is known to account for the single particle distributions in e + e~ collisions, 
we fit the rapidity and transverse mass spectra of pions and kaons measured in Pb+Pb collisions 
by NA49 collaboration. This model nicely fits both the rapidity and the transverse mass spectra of 
each particle species with a few parameters. However, fitting all the particles with a single value of 
T p is not possible. This analysis shows that the success of the thermal models does not necessarily 
mean the thermalization of the system in collision. 



PACS numbers: PACS numbers: 97.60. Jd,21.65.+f,12. 38. Mh 
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I. INTRODUCTION 

It has been of great interest to find out whether a 
hot and dense nuclear matter has been formed during 
the high energy nuclear collisions 1] , and when the tem- 
perature and/or density is high enough one expects a 
system of free quarks and gluons, namely quark-gluon 
plasma. Main issue here is the thermalization of the par- 
tonic or hadronic system formed just after the collisions 
and throughout the evolution of the following hadronic 
system. Thus probing the thermalization of the system 
in collision is of great importance. The large elliptic flow 
and the jet quenching effect indicate a collective behav- 
ior. However those are not yet the direct experimental 
evidence of thermalization. Best one can do so far is 
to assume thermalization and to show the experimental 
data is not inconsistent with the assumption of thermal- 
ization. In this case one has to show that assuming ther- 
malization is the only way to explain the data in order to 
claim the thermalization. We will show in this paper that 
a simple statistical model can fit both the rapidity and 
the transverse momentum spectra of various hadrons. 

In the relativistic heavy-ion collisions at AGS, SPS 
and RHIC, hadrons are measured quite accurately, espe- 
cially for pions, kaons and protons. Usually the hadronic 
data are presented in many different ways: (1) total 
multiplicities of particles or the ratios of the multiplic- 
ities of different hadrons, (2) single differentials such 
as the rapidity distribution dN/dy, and the transverse 
mass spectra dN/m±dm± , (3) the double differentials 
such as B-E correlation function and (4) the averages 
or moments of certain physical quantities such as fluc- 
tuations and flow coefficients which is the second mo- 
ment of the azimuthal angle, 4>. The total number of 
particles can be obtained either by integrating the rapid- 
ity distribution over the whole range of the rapidity, i.e. 
Ni = J dy(dNi/dy) or by integrating over the whole p± 
range, i.e. Ni = J m±dm±dN/m±dm±. 

Models assuming thermal and chemical equilibrium fit 
each pieces of data very nicely 0, El IE IE an d ex- 
plain many features of data. Thermal models assuming 
chemical equilibrium 0, 0, 3 , when applied to the ratios 
of various hadron multiplicities fit the data impressively 



well with only a few parameters such as the temperature 
T c , baryonic chemical potential, fib and/or the strange 
chemical potential, fj, s . 

Fireball models with strong radial flowJEIEQ] fits the 
transverse mass spectra quite nicely with a few parame- 
ters, Tth, fJ-b and radial expansion velocity at the surface, 
ft. 

The rapidity distributions are usually Gaussian. 
Spherically expanding fireball model applied to pion data 
at AGS energyQ results in too narrow width compared 
to the data, and in the cylindrically expanding fireball 
model|E| fast longitudinal expansion makes the rapidity 
width wider and seems to fit the data when the system 
sizes are small. However, expanding fireball models with 
any geometry seems to have difficulty in fitting the ra- 
pidity distributions of pions and kaons, simultaneously, 
in the Pb+Pb collision at SPS. 

The two different sets of freeze-out parameters from 
the ratio analysis and the transverse momentum spectra 
analysis leads to the conjecture that the chemical com- 
ponents freezes earlier at high temperature, T c h, and the 
thermal equilibrium is still maintained until T t h, keeping 
the number of certain particle species constant 8]. 

However, calculation by Broniowski and FlorkowskiQ 
claim that both the ratios and transverse mass spectra 
can be fitted with a single set of freeze-out parameters. 
As the models seem to be the similar, independent cal- 
culations are needed to check which one is right. 

In spite of the success of thermal models for hadrons 
produced in the relativistic heavy-ion collisions, it is 
known that the statistical model of m ultip article pro- 
duction leads to the Boltzmann factor |9L fiot tll| and can 
explain the hadronic data produced in e + e _ and pp col- 
lisions. In the series of papers by T. T. Chou, C. N. Yang 
and E. Yen a statistical model for multiparticle pro- 
duction is derived and applied successfully to the rapid- 
ity distribution and transverse momentum distribution of 
measured hadrons in c+ c~ collisions by TASSO collabo- 
ration. This model has been further modified and studied 
by Hoang [l^ to include the Lorentz-boosted Boltzmann 
factor in order to mimic the shift of the peak of the ra- 
pidity distribution in the asymmetric collisions. It should 
be emphasized that this statistical model does neither re- 
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quire the thermalization of the system nor the concept 
of freeze-out. Recently the same argument has been re- 
vived 0,0,0], but application of the statistical model 
to the relativistic heavy-ion collision data has not been 
attempted yet. 

Thus it is very interesting to study whether the sta- 
tistical model which successfully fit e+e - collision data 
can also fit the rapidity distributions and transverse mo- 
mentum spectra of pions and kaons produced in the rela- 
tivistic heavy-ion collisions, especially whether it can fit 
all of them with a single set of parameters. A priori it 
is not clear whether pions and kaons have the same pa- 
rameter "T" in the statistical model. However, there is 
a claim that e+ e~ collision data could be fit with only 
one "temperature" for various hadrons [ill 

In this paper we will show that a simple statistical 
model can describe both the rapidity and transverse mass 
spectra of various hadrons measured in the heavy-ion col- 
lisions. In our previous work|18j. the rapidity distribu- 
tion of hadrons are analyzed using the same model. 

In the next section the statistical model by chou, Yang 
and Yen is reviewed and in sec. Ill the result of our 
analysis of the rapidity distributions and transverse mass 
spectra of 7r + , 7r~, K + and K~ in Pb+Pb collisions at 158 
GeV- A by the NA49 collaboration^^ E3 is presented 
and finally we summarize in sec. IV. 



II. STATISTICAL MODEL 

In this section we follow the paper by Chou, Yang and 
Yen ^l| for the statistical model invented for e + e~ colli- 
sions which can be used for the heavy-ion collisions with 
minor modifications. 

We will assume the multiparticle production in the rel- 
ativistic heavy-ion collisions is stochastic but subject to a 
number of conditions: (a) energy conservation, (b) lead- 
ing baryon effect which means that leading baryons in 
both the projectile and target region do not participate 
in the collision, (c) d 3 p/E probability for each particle, 
and (d) transverse momentum cutoff factor g(p±). Then 
the probability for non-leading particles in the collision 
zone can be taken as 



5 ( E i ~ E h ) J] (d^/E,) g{p u ) 



(1) 



where Eq = i/s and h is the fraction of the total energy 
Eq used for the particle production in the central region. 

It was pointed out that the probability distribution is 
exactly the same as in the microcanonical ensemble in 
the statistical mechanics. By the well-known Darwin- 
Fowler method the single particle distribution of such 
an ensemble can be converted into a distribution in the 
canonical ensemble : 



Probability = (d 3 p / E) g(p ± ) e^p(-E/T p ) (2) 



where T p is called the partition temperature. Noting that 
d 3 p = dpLp±dp±d<fi = Edyp±dp±d<p, one gets the equa- 
tion for the rapidity distribution by integrating over the 
transverse momentum. 



dN l /dy = 2irK p±dp± gi(p±)exp(-Ei/T p ) (3) 

Jo 

where if is a normalization constant, and p max = E h 
which can be taken as infinity in the relativistic heavy- 
ion collisions. E{ — nir cosh y , and my = yj mf + p\ . 
We will use for the high momentum cutoff factor as in 
Ref. [ill, 



9i(p±) = exp(-aif>jj 



(4) 



where a* related to the average transverse momentum, 
< p± >. One can show that for the massless case, 



at = 2[<p u >] 1 



(5) 



and in fitting e + e _ data, the value of a was taken from 
the measured data. 

Similarly the transverse momentum or transverse mass 
spectra is given, after integration over the rapidity y, as 

dNi/m±dm± = 2TrKgi(p±) J dyexp(-Ei/T p ) (6) 

Eq. © can be integrated for the massless particles and 
the role of parameters can be easily understood. 



dNJdy 



2ttK 



(a, +coshy /Tp) 2 
whose maximum at y = is 

2irK 



(7) 



(dNi/dy) r 



2nKT 2 

(8) 



{a, + 1/T p ) 2 { ai T p + 1) 



and the width Syi at the half maximum is given from the 
relation 



coshSyt = (V2 - l)a, t T p + V2 



(9) 



For the massless particle the width is governed by the 
factor a^p. Thus when the value a is taken from ex- 
periment, the parameter T p is a measure of the width 
of dN/dy in this model. This is quite different from 
the thermal models, where T p determines the slope of 
the transverse momentum spectra and in order to fit the 
width of the rapidity distribution one further need the 
concept of longitudinal expansion of the system. 

This can be understood again from Eq. @, which de- 
scribes an exponential shape. Contrary to the naive ex- 
pectation, the slope of the exponential is determined not 
by T p , but by 1/a through Eq. (@J. This is because the 
fitted value of T p is very large ( ~ 1 GeV) compared to 
the average transverse momentum, which usually trans- 
lates into the temperature. Thus in this model the main 
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slope of the exponential shape is put in by hand as a high 
momentum cut-off factor, Eq. Q and the integration in 
Eq. @ gives a modification to the simple exponential 
shape. 

In Eq. Q and Eq. Q , K, T p and oti 's are parame- 
ters. In the next section results of fitting of the rapidity 
distributions and transverse mass spectra of 7r+, tt~ , K+ 
and K+ in Pb+Pb collisions at 158 GeV- A by the NA49 
collaborations. 



III. ANALYSIS 

Main results of the fitting of the rapidity distributions 
and transverse mass spectra of 7r + , tt~, K + and K _ 
in Pb+Pb collisions at 158 GeV-A measured by NA49 
collaboration^! 0] are shown in Fig. 1-4 and the re- 
sulting parameters are tabulated in Table 1. 

In Table 1, we have fitted each species separately and 
thus the resulting parameters arc all different. < p±_ > 
is calculated from a using Eq. (j5j). One may use exper- 
imental value of a and then K and T p are the only pa- 
rameters in this model. It is amazing that a very simple 
model with only two parameters can fit both of the ra- 
pidity and transverse momentum spectra of each particle 
species. 

In this model there is no a priori reason for the com- 
mon value of T p for different particle species. Simulta- 
neous fitting of the rapidity and transverse mass spectra 
of pions and kaons with a single value T p is not possi- 
ble, which is contrary to the case of e + e~ collisions [l2|. 
However, in Ref. |l2( Lorentz-boosted exponential was 
used and direct comparison is not possible and further 
studies with the Lorentz-boosted exponential is needed 
to check whether in the relativistic heavy-ion collisions a 
single value of T p can describe the particle spectra of all 
the different species. 

In Ref. 0| we have successfully fitted only the rapidity 
distribution of pions and kaons with a single T p but differ- 
ent ati's. We get K = 1212, T p = 0.98 GeV, a w + = 5.98 
c/GeV, a n - = 5.67 c/GeV, a K + = 11.4 c/GeV, and 
a K - = 17.1 c/GeV. The large values of oti's spoils the 
interpretation of a in terms of the average transverse mo- 
mentum through Eq. (jSJ). This means that even though 
all the rapidity distributions are fitted simultaneously, 
the transverse mass spectra are far off. However, one can 
also fit the rapidity distributions of each particle species 
separately by fixing the cti values as the experimental val- 
ues and the resulting values from the least square fit are 
approximately the same as those in the Table 1, meaning 
that by fitting only the rapidity distribution one can fit 
the transverse mass spectrum also. Concentrating only 
for the rapidity distributions, one can either fit them with 
a single value of T p with the transverse mass spectra far 



off or fit both of the rapidity and transverse mass spec- 
trum with different T p values for different particle species. 

As was discussed in the previous section T p is above 1 
GeV and thus in Eq. © , 1/a plays the role of the slope 

TABLE I: The fitting parameters of the rapidity and trans- 
verse mass spectra 





K 


T p 
GeV 


Q 

c/GeV 


<P±> 
GeV/c 


X 2 /n 


7T+ 


599.9 


1.27 


4.23 


0.47 


5.4 


7T _ 


649.8 


1.19 


4.15 


0.48 


15.0 


K+ 


66.86 


1.57 


2.69 


0.74 


2.5 


K' 


37.82 


1.35 


2.79 


0.72 


2.7 



parameter and the integral in Eq. © acts as the high 
energy cut-off. 

For pions we have used K, T p , a n + and a n - as parame- 
ters. One may use the same T p and a but different K val- 
ues for 7r + , 7r~. Or one may further introduce the charge 
chemical potential fi c to write K v + = K exp (fi c /T p ) and 
K^- = AT exp (—fj, c /T p ). In this way one gets K = 625.0, 
T p = 1.21 GeV, a = 4.17 c/GeV and ^ c = 0.043 GeV. 
The quality of the fitting is more or less the same. Sim- 
ilarly for kaons one can introduce the strange chemi- 
cal potential fj, s to write K K + — K exp (— /i s /T p ) and 
Kk- = K exp (fj, s /T p ). Then the rapidity and transverse 
spectra of both the K + and K~ can be fitted simultane- 
ously with K = 50.7, Tp = 1.43 GeV , a K+ = 2.6c/GeV, 
a K - = 2.9 c/GeV and (i s = -0.415 GeV. 



IV. CONCLUSION 

We have fitted the rapidity and transverse mass spec- 
tra of 7r + , 7r~, K + and K~ in Pb+Pb collisions at 158 
GeV-A measured by NA49 collaboration with a simple 
statistical model. For each particle species both the ra- 
pidity and the transverse mass spectra are nicely fitted 
with only three parameters. However, simultaneous fit- 
ting for all the pions and kaons with a single value of T p 
is not possible. 

Thus the success of thermal models for the hadron mul- 
tiplicities and/or transverse momentum spectra, by itself 
alone, does not necessarily mean the thermalization of 
the system. 
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FIG. 1: Rapidity distribution and transverse mass spectra of 
7T + in Pb+Pb collisions at 158 GeV-A by NA49 collaborations. 
Lines are the fittings with Eq. and Eq. © with parameters 
in Tab. (1). In the right figure from the top y bin is 0.0 ±0.1 
and the second one is for y =0.2± 0.1, and so forth. 
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FIG. 4: Same as in Figl. for K" 



